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Abstract 

Based  on  the  Navier-Stokes  equations,  the  compressible  viscous  vortex  equation  des¬ 
cribing  the  separated  vortex  of  an  aircraft  is  derived  and  a  numerical  method  for  sol¬ 
ving  this  equation  is  developed.  It  can  be  used  to  subsonic,  transonic  and  supersonic 
vortex  flown  The  breakdown  condition  of  the  oompressible  separated  vortex  is  derived 
and  used  to  compute  the  breakdown  place.  A  series  of  interesting  features  about  the  ef¬ 
fects  of  various  parameters  such  as  compressibility  etc.  on  separated  vortex  motion  arc 
obtained  and  analysed. 


AN  INVESTIGATION  ON  THE  COMPRESSIBLE  VISCOUS  SEPARATED  VORTEX  OVER 
THE  WING  LEADING  EDGE 


Lin  Bingqiu 
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< 


Based  on  the  Navier-Stokes  equations,  the 
compressible  viscous  vortex  equation  for  de¬ 
scribing  the  separated  vortex  of  an  aircraft 
is  derived  and  a  numerical  method  suitable  for 
subsonic,  transonic  and  supersonic  flows  is 
given  which  can  compute  the  characteristics  of 
separated  vortex  or  trailing  vortex.  The  break¬ 
down  condition  of  the  compressible  separated 
vortex  is  derived  and  used  to  compute  the  break¬ 
down  point/  A  series  of  interesting  features 
about  the  effects  of  various  parameters  such  as 
compressibility,  etc.  on  separated  vortex  motion 
are  obtained  and  analyzed. 


I .  Preface 

The  incident  angle  adopted  by  recent  winged  aircraft  is 
becoming  larcer  and  larger  in  order  to  more  effectively  utilize 
the  vortex  lift  generated  by  separated  vortex  of  the  component. 

But  along  with  it  comes  the  interference  between  the  vortex  system, 
between  separated  vortex  and  component,  and  between  aircraft 
and  trailing  vortex,  and  the  future  motion  state  of  the  vortex 
itself  also  has  unpredictable  factors  occurring  sometimes.  Therefor 


a  detailed  study  of  the  motion  laws  of  the  vortex  itself  becomes 


the  first  step.  Up  until  now,  there  has  been  much  research,  both 

theoretical  and  experimental,  done  on  this  low  speed,  incompressible 

Til  r  2 1  r  3 1 

vortex.  Brown  ,  Hall  ,  and  Muirhead  ,  etc.  have  done  much 
research  on  the  compressible  viscous  vortex,  but  it  is  still 
far  from  sufficient. 

The  effects  of  compressibility  on  the  vortex  are  extremely 
distinct.  We  must  establish  a  computation  method  suitable  for 
subsonic,  transonic  and  supersonic  flows  and  clearly  understand 
the  effects  of  various  parameters  as  well  as  the  breakdown  features 
of  the  compressible  vortex,  etc.  This  is  exactly  the  subject 
which  this  paper  wants  to  solve. 

II.  Motion  Equation  of  the  Separated  Vortex 

We  assume  the  vortex  is  laminar  and  symmetrical  and  can 
be  described  by  N-S  equations  using  cylindrical  coordinates. 

The  coordinates  system  x,  r,  9  correspond  to  component  velocities 
u ,  v ,  w .  The  external  flow  (or  external  flow  field)  of  the  vortex 
can  be  described  by  compressible,  nonviscous,  isentropic  flow. 
Assuming  the  vortex  is  very  slender,  let  £  be  the  local  slender 
ratio  of  vortex,  and  £  is  first-order  small  quantity.  Based 
on  quantity-net  analysis,  v=  £ u  is  derived  from  the  continuity 
equation;  we  take  w=0(u).  The  coordinates  and  various  physical 
quantities  are  written  as  normalized  quantity-net  forms  according 
to  this  order  of  magnitude  relations:  r=r/R^,  x-£x/R^,  u=u/u^, 
w  —  w  / 1  j  ^  ,  v  =  v/£  u.,  p  =  p/  ^  ^  u?  ,  p=p/p  ,  T  =  c  T/uT  and  h  =  h/uT  . 
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The  bar  above  the  nomenclatures  denotes  normalized  quantity- 


net  quantity.  R  ^ ,  p  ^  aid  u^  are  initial  radius,  density  and 

axial  velocity  at  the  outer  edge  of  the  vortex  cross  section, 

respectively,  and  here  they  are  selected  as  reference  quantities. 

p,  T,  h  and  c  are  pressure,  temperature,  enthalpy,  and  specific 
P 

heat  at  constant  pressure,  respectively.  Substitute  the  above 
normalized  quantity-net  quantities  into  the  N-S  equations,  and 
take  the  reciprocal  of  Reynolds  number  R(  3=  P  iUiV/^  ^  as  a  first- 
order  small  quantity,  then  omit  the  second-order  small  quantities 
in  the  equations  to  obtain  vortex  motion  equations  that  describe 
the  compressible  viscous  vortex: 

continuity:  d(r<??2.  +  Mlfi- 1  —  p  (1) 

dx  dr 


Axial  Momentum:  au  +  pv  —  —  _  h - L_ +  _L  JL  (  / 

dx  Or  dx  rRc  dr  Re  dr  V  dr  J 

Radial  Momentum:  t,u,i 

.r  dr 

»  '  ‘  -  •  \  •- 

Rotational  Momentum:  ^  dw  t  fll.  dw  (  pwv  __  2  f  /  dw  «y\ 

dx  dr  r  .  rRe  'dr  r  r 

■Re  dr  l  A  Sr  r  /  j 

Energy  Equation:  „  Q  +  dT  _  u  d£  +  v  df_  +  J_\(  d»Y 

dx  dr  dx  "dr  >  l 'Or  / 
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(4) 


(  5  ) 


where  all  quantities  are  normalized  quantity-net  quantities, 
and  the  bars  above  nomenclatures  are  omitted  for  the  sake  of 
convenient  writing.  is  the  Prantdl  number,  fA  / JA,  .  ,  the 

equation  groups  is  a  parabolic  type ,  and  solutions  can  bo  derived 
if  initial  and  boundary  values  are  known. 


III.  Numerical  Method  of  Vortex  Equation 


The  method  for  solving  partial  differential  equations  (1)~ 

(5)  should  be  based  on  the  mission  of  the  problem.  The  net-grid 
method  is  suitable  for  computing  vortex  features  in  a  certain 
section  of  the  flow  field,  such  as  computing  the  features  in 
the  breakdown  region.  Even  so,  convergence  is  difficult  to  obtain 
at  a  slightly  higher  Reynolds  number.  It  is  not  economical  using 
this  method  to  compute  a  slender  trailing  vortex  either.  In 
order  for  convenient  theoretical  investigation,  we  turned  to 
the  solution  method  similar  to  those  adopted  by  Reference  [4] 
and  [5],  and  expanded  it  to  solve  for  vortex  motion  taking  compress¬ 
ibility  into  account. 


1.  Solve  for  the  radial  velocity  v  from  equation  (1), 

„  _  _  _L  [  Kisul  dr, 

rp  J»  dx 

and  substitute  into  the  v's  in  equations  (2),  (4),  (5).  Then 
take  the  integral  along  the  radial  direction  while  noting  the 
outer  edge  and  core  condition  of  the  vortex 

du  i  du  |  n 

Differential  and  integral  equations  corresponding  to  (2),  (4), 

(5)  can  be  obtained: 

(  ~~  l  rPu  (“  —  vr)]dr  +  l"  rpudr  —  —  (  r  4r  (  6  ) 

Ox  dx  Jo  !«  dx 

(  [  rpu(c  —  C)]  dr  +  —  f  rpudr  —  —  Sit  (-) 

dx  dx  Jo  Re  1  '  ' 

(  [rpu  ( T  —  Tr)]dr  +  (  rpudr 

Jo  dx  dx  J« 
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where  c-rw  and  C  =  Rwe  are  vortex  internal  and  outer  edge  circulation 
parameters,  respectively,  and  represent  vortex  intensity.  The 
subscript  "e"  represents  the  vortex  outer  edge. 


2.  Expand  the  functions  u,  w,  p  inside  the  integral  sign 
into  polynomials: 

m  «  * 

“  —  s  “<v.  * — yi  w>' i%  p  ~  s 
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where  =r/R  and  the  coefficients  in  the  expansion  are  determined 

by  boundary  conditions.  If  u,  w,  p  are  known,  then  p  is  obtained 

.IT  P 

from  (3)  and  T-  is  obtained  from  the  state  equation  T=  ^  p  in 
normalized  quantity-net  form.  < 


At  the  outer  edge  of  the  vortex  where  r=R,  u=u  and  — 

<=  dr 


d>u 

5r‘ 


CIR, 


and  ^  and  &  „ 

dr  dr  Qr  R,  • 


S  w 

At  the  vortex  core  where  r  =  0,  u  =  u  ,  p  =  p  ,  w=0  and  —z — =w. 

o  o  or 

w  is  called  the  vorticity  parameter  of  core.  Its  relations  with 

the  vorticitv  of  the  vortex  core  is  0  =2w/R;  and  also  from  symmetrical 

o 

feature  dp  _  du  ^  dp  _  0 

dr  dr  dr 


where  — <j~r  '  J  r  _q  =  0  is  an  additional  condition  for  p.  Each  expansion 
has  5  boundary  conditions;  therefore,  pentanomials  are  obtained 
as  follows: 


«  —  «,(.!  —  6V  8  V  —  3  V)  4-  »,(6n>  —  *V  +•  3 VI 

if  —  —  (10—  1 5ij  +  6V)V  +■'  1  —  Ji)  +  3*|*  — 

R 
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r_ 
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The  subscript  "0"  denotes  core  vorticity. 


3.  In  principal,  substituting  (9)  into  each  differential 
and  integral  equation  ( 6 ) — ( 8 )  and  then  taking  the  integral  sequential 
can  simplify  them  into  a  group  of  differential  equations.  However, 
since  the  form  of  integral  function  is  more  complicated,  it  is 
extremely  cumbersome  to  derive  using  an  analytical  method.  For 
this  reason,  we  adopted  a  more  generalized  method.  Separate 
the  first-order  derivatives  from  the  equations,  then  set  the 
integral  of  various  complicated  forms  as  the  coefficients  of 
those  first-order  derivatives,  and  use  a  numerical  method  to 
compute  these  coefficients.  In  order  to  illustrate  the  derivation 
method,  a  term  with  the  most  representative  form  is  selected 
for  demonstration  purposes,  for  instance: 


where 


separate 


<*  —  <£(>,(*),  *•(*),  R(*),  «>(x), />,(*)»  “,(*X  COO) 

^  —  'KPoC*),  «.(*),  R(*),  w(x),  />,(*),  «,(*),  c 00) 

ail  the  derivatives,  then  the  above  integral 


is  written 


a  s 


f  “  (S  (iy  if  y  5  f'y; 


W=u  ,  Y„  =  C. 

o  C 


In  n  and  m<0, 


Y  =n  ,  Y_=u  ,  Y o  =  w ,  Y  .  =  R  ,  Y_  =D  , 

1  ^  o  2  o  3  4  d  •  c  ul. 

an  odd  point  will  appear  at  the  lower  limit  of  the  integral  for 
the  function  under  integration.  But  the  odd  point  can  be  eliminated 
within  the  range  of  the  problems  we  are  solving. 


Now  three  first-order,  near-linear  differential  equations 

are  obtained,  vet  there  are  four  functions  p  ,  u  ,  R,  w  to  be 

o  o 

solved.  Use  a  method  similar  to  thos  1  in  Reference  [4]  and  [5] 
to  make  an  equation  (2)  satisfy  the  boundary  condition  at  the 
core  and  obtain  a  supplementary  equation.  Thus,  a  closed  equation 
system  is  formed: 


-fft  V 

«„ 

R 

L  u>  -» 


-  B 


r  p. 

i 

c 

L  0  J 


■VT" 
U 
U 

■U- 1 


( 10) 


Ihe  apostrophe  means  to  take  the  derivative  once  with  respect 


to  x,  where 


au  ■+■  a  „  +  av  an  +  tf>« 

bu  b,  j  bn  b  It 

c„  -  t.,  +  c.i  fu  +  r,j  c„  +  c.i  —  Cy,  4-  c,?  '  c „  —  c -h  r,. 
Jn  d*  0  0 

r —  (<*,2  +  o.:)  —  (fli.  +  «ji  ~t-  o»)  +  a»  0 

—  bu  —bu  —(b„  +  b3l)  0 

—  cu  —  —  (r„  +  ca  —  +  <■«)  —  Ci*  —  —  0 

0  0  0  0 


See  appendix  for  the  expressions  of  the  aforementioned  elements. 

- .  Equation  (10)  is  a  differential  equation  with  functions 
of  o  ,  u  ,  ?. ,  w  and  with  the  x-axis  coordinate  as  so  1  f  -  va  r  i  a  nee . 


the  vortex,  axial  velocity  u  ,  circulation  parameter  C,  Reynolds 

number  Rg,  Prantdl  number  P  and  viscosity  coefficient  ratio  f. 

It  reflects  physically  the  effects  of  nonuniform  external  flow, 

viscous  diffusion  and  heat  conduction  on  vortex  motion.  The 

initial  conditions  at  the  outer  edge  of  the  vortex  are  uc  =1  and 

i 

R^  =  l."  The  initial  pressure  of  external  flow  can  also  be  expressed 

in  the  form  of  external  flow  initial  Mach  number  M  .  M  =[(1-C“)/ 

a .  a .  1 

i  i 

1/9 

/  pc  ]  “  .  In  the  computations  hereinafter,  P  and  f  are  approxi- 

1  n  7  s 

mately  taken  as  Pr=l  and  f=(T/T^)  *  .  When  the  initial  values  of 

external  flow  parameters  and  functions  are  given,  then  the  Rung-Kutta 
method  of  fourth-order  accuracy  can  be  used  to  numerically  solve 
equation  (10).  It  is  not  limited  by  external  flow  M  number. 

IV.  Analyses  and  Discussions  of  Calculated  Results 

1.  Typical  Internal  Distribution  Features 

W'e  selected  a  group  of  initial  parameters  and  external  flow 
parameters  which  correspond  to  stable  motion  of  separated  vortex,  and 
the  calculated  results  are  as  shown  in  Fig.  1.  It  displays  the  radia 
changing  laws  of  various  physical  quantities  on  two  cross  sections  of 
the  vortex  in  supersonic  external  flow. 
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Effects  of  compress ibilitv  on  initial  pressure  at  the  vortex 


In  the  initial  stage,  x=1.65,  and  the  changes  in  various  physical 
quantities  along  the  radial  direction  are  greater.  Further  downstream, 
they  then  tend  to  stabilize  and  should  become  equal  to  those  of  the 
external  flow.  For  every  cross  section,  axial  velocity  u  near  the 
v o  r  d  c ora  i  3  ma  ima  1 ,  f.;h o  r  o  a  3  d  o n  3  i  d  v ,  n  r  a  .=  3  *  1  r  ^  m  ;i  t  or>  o o  r  V:  *  i  r  ■-*  I  r 
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the  opposite.  The  radial  velocity  v  basically  maintains  its  value  as 
a  second-order  small  quantity,  which  is  consistent  with  the  order  of 
magnitude  analysis.  At  the  vortex  core,  axial  velocity  is  in  the 
deceleration  process  and  gradually  approaches  1.  Pressure  change 
is  in  the  inversed  pressure  process.  Rotational  velocity  w  is  similar 
to  Berger  distribution;  it  is  0  at  the  vortex  core  and  after  reaching 
maximum  (  ^  <  1)  as  ^  increases,  it  approaches  0  at  infinity.  The 
vortex  radius  gradually  expands  and  this  is  closely  related  to  vortex 
core  deceleration. 

2.  Compressibility  Effect 

# 

The  internal  density  of  the  vortex  under  high  speed  motion  is 
variable.  But  p  <.  0  must  be  maintained;  it  is  the  necessary  physical 
condition  for  the  vortex  to  be  in  existence.  Substitute  the  expansion 
of  p  into  (3)  and  let  r=0,  then  the  initial  density  at  vortex  core  is: 


fherefore,  if  satisfied 


p„  -  po,  >  cm 


(in 


or 


Po,  < 


1  +  C]  c] 


TMtfj  A 

Thus  ,  Po.>  0  can  be  guaranteed.  Equation  (11)  shows  that  only  if  a 
u  i 

specific  pressure  difference  is  guaranteed  can  be  existence  of  the 
vortex  be  maintained,  and  that  at  this  time,  the  pressure  at  the  vortex 
core  must  be  lower  than  the  external  pressure.  Equation  (12)  shows 
that  the  selected  initial  pressure  and  vortex  intensity  for  the  exis¬ 
tence  of  the  vortex  are  related  to  external  flow  Mach  number.  This 
relation  is  visually  displayed  bv  Fig.  2.  Solid  lines  are  the  p 


curve  at  constant 


curve  at  p  =0  and  constant  C..  Dotted  lines  are  p 

*0 .  1  c . 

l  i 

C^.  The  region  between  solid  lines  and  straight  line  pQ  =0  (the 

abscissia)  is  the  region  where  significant  vortex  core  pressure  p  is 

u  i 

selected,  i.e.,  the  region  where  the  vortex  can  be  formed.  As  Ma^ 

increases,  said  region  becomes  smaller.  After  Ma^  is  greater  than  1, 

this  effect  is  more  distinct.  It  stabilizes  after  Ma .  >  1 . 7 .  The 

i 

larger  Ch  is,  the  smaller  the  region  is.  We  have  noted  that  when  is 
larger,  and  after  Ma ^  exceeds  a  certain  limit,  said  region  shrinks  to 
zero.  At  this  time,  it  is  impossible  to  form  the  vortex.  Moreover, 
under  the  condition  of  smaller  CM  (as  CM =0.1  in  the  figure)  and  very 
large  Ma^,  the  pressure  difference  between  the  inside  and  outside  of 
the  vortex  is  almost  negligible.  At  this  time,  it  is  meaningless  to 
generate  lift  using  separated  vortex.  In  experiments,  the  measure¬ 
ments  of  pressure  distribution  for  a  slender  cone-shaped  lift  body 

r  6 1 

have  shown  1  that  the  pressure  peak  on  the  surface  trailing  the 
wind  becomes  level  as  the  supersonic  Ma  number  increases,  which  means 
the  disappearance  of  the  effects  of  separated  vortex. 

3.  Breakdown  Phenomenon  of  Compressible  Vortex 

We  assumed  that  there  is  a  stationary  point  at  the  vertex  core, 
and  the  origin  of  coordinates  was  set  at  the  stationary  point  with 
positive  direction  of  x-axis  pointing  downstream.  Each  physical 

quantity  is  expanded  near  the  x-axis  of  the  stationary  point: 

"  -  «,<»  4-  «,OV  +  ►}(*) r3  +  ... 

u  —  w,(x )r  4-  K'jOV  4-  •  p  —  f„(*)  4-  p,OV  4-  •  •  •  > 

P  “  Pt 0*0  4-  4-  •  ■  • 

Substituting  info  the  X-S  equations  and  take  ing  1/Pe  as  higher  than 
second-order  small  quantity,  then  the  following  is  obtained: 


; 
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p,«0  +  2p0t',  —  0 


Pt  +  P&<>u't  —  0 


uaw,  +  2»<1«p1  —  0 


From  (14)  and  (15)  to  obtain,  respectively 


t-  -  _  1  • 

t>,  — - u 


Substitute  (17)  into  (16)  to  obtain 


u,  —  e,iv, 


where  c^  is  an  integral  constant.  That  axial  velocity  should  have  the 
following  features  near  the  stationary  point  are  already  known: 

“i  r-c~  <  0  or  (“5)'U-«- <  o,  u;|r_+ <  0  or  («D'lx-,*>0 
Therefore,  it  can  be  observed  from  (18)  that  pQ  increases  from  up¬ 
stream  of  the  stationary  point  to  the  stationary  point  and  pressure 
is  inversed  along  the  axis.  Conversely,  pressure  changes  along 
the  stream  from  the  stationary  point  to  downstream,  i.e.,  pressure 
is  maximal  at  the  stationary  point.  It  can  be  derived  from  (19) 
that  w1  is  egual  to  zero  at  the  stationary  point  and  that  in  the 
upstream  and  downstream  areas  near  the  stationary  point,  w^  changes 
sign  after  passing  the  stationary  point  because  uq  changes  sign. 

This  shows  that  at  the  stationary  point  the  radial  gradient  of 
rotational  velocity,  ^w/  r  I  r_0~0 ,  or  vortecity  Q  =  0  (or  w=0). 

After  the  vortex  passes  the  stationary  point,  the  rotational  velocity 
of  vortex  reverses  direction.  Equation  (17)  shows  that  the  radial 

gradient  of  radial  velocity,  »  _  )>  0 .  Therefore,  there  is  a 

c?  r  1  r  —  u 


neqa  1 1 vo 


;h  surface  stating  from  the  stationary  ooint.  The 


flow  departs  from  the  vortex  core  with  the  stationary  point  as 

the  separation  point,  and  it  is  also  the  odd  point.  Its  mathematical 

condition  is  concluded  as: 

-0,  f4  >0  (20) 

dr  r*  dr 

It  is  the  axially  symmetrical  breakdown  condition  of  the  vortex. 

It  is . completely  identical  in  form  to  that  of  the  incompressible 

„  [4] 

vortex 

A  large  amount  of  computation  and  experiments  for  subsonic 
and  transonic  vortex  flow  have  shown  that  the  breakdown  point 
takes  two  forms  due  to  the  effects  of  compressibility.  One  is 
the  sudden  change.  Immediately  upstream  of  the  breakdown  point, 
various  physical  quantities  gradually  change.  At  the  breakdown 
point,  the  realization  condition  (20)  can  be  approached.  Figure  3 
shows  the  calculated  results  under  uniform  external  flow  and  subsonic 
condition.  Breakdown  place  XBD=10.2,  and  uq  and  w  approach  zero. 

We  studied  the  effects  of  various  parameters  such  as  Ch  , 
u  .,  Re  and  initial  pressure  difference  /i  p.  =p  .  -  p  .  on  the  break- 

down  point  of  the  vortex  through  a  large  amount  of  computations. 
Increasing  Ch  ,  Ap.  ,  respectively  or  reducing  uq^,  Re,  respectively, 
will  promote  breakdown  of  the  vortex.  Conversely,  they  will  delay 
the  occurrence  of  breakdown.  These  conclusions  are  completely 

r  4  l 

identical  to  the  results  of  the  incompressible  vorte.x 


Fig.  3.  Computation  of  the  breakdown  point  of  subsonic  separated 
vortex 


■\  c  p  e  n  d  i  x 


a, i  —  R‘  7  Jj-  C>  —  1»  2,—,  6) 

J|t  =  R'J'  7  -jL-  Cp«(a  -  a,)V7  +  2rJ ^7|0“(«  "  “«V7 

J„  =  R‘  J4  7  pudqx  au  =  R'  7  <<7  (■  =  l>2,  4,5,'5) 

=  R*I1  7-|f-  drJJrlR  j'^n-RP. 


k^R‘\\r,  d^-£l  Jr)  <i-l,2,.-,6) 


dpuQc  —  C) 


a« 


^7  +  2R  ||  7Ru(f  -  C~)Jn,  b„  =  R'  js  W7 


<•„  -  3. 5R‘  ||  7-^-  C-O  -  0.25T.p)y7  +  2R  j#  7-(P  -  0 .  l5T.p)J  7 

(i  -  l,  2,  -,  S) 

f „  =■  —  f  r<b/7»  fu  *»  R'  f  7“^7  (•=•!»  2,<,  I,  6) 

e-  - *'\l.  * 4*  dr> -* 

e,t  mm  R*  ff-^~  P.<^7  2,  —  ,7),  rf.i  I»  dll  ”  ■ 

».  -  (-£)■.  ^4r  1:  '(£ ■© . 
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Abstract 
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'  j^The  principle  of  flow-yelocity  ,xnd  flow-pattern  measurement*  using  Doppler- frequ¬ 
ency-shifted  absorption  (or  emission)  spectral  lines  as  well  as  Doppler-frequt^  y-shif- 
ted  .and  deformations]  non^saturation  spectral  line-shapes  is  discussed,  ia— this  paper. 
For  the  case  of  gas  flow  with  simultaneously  inbomogeneausly  and  homogeneously  bro¬ 
adened  line,  the  general  expressions  of  non-saturation  spectral  line-shape  and  Doppler- 
frequency -shifted  spectral  line  are  presented.  For  the  case  of  three  representative  gas 
flow  with  inhomogeneously  broadened  line,  thd  oancrete  expressions  of  non-saturation 
spectral  line-shapes  and  the  qualitative  and  quantitative  relationships  among  flow-velo- 
eity,  flow-pattern,  Doppler -frequency-shifted  spectral  line  and  defonnational  non-satu- 
ratioa  spectral  line-shape  are  also  given. 
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This  paper  investigates  the  problem  of 
using  frequency-shifted  and  deformational 
laser  spectral  lines  to  measure  the  gas  macro¬ 
scopic  flow-velocity  and  f low-pattern .  For  the 
case  of  simultaneous  effect  by  broadened  Doppler 
lines  and  broadened  pressure,  the  general  expres¬ 
sions  of  nonsaturation  spectral  line-shapes  and 
frequency-shifted  relationship  are  obtained;  for 
the  case  of  three  typical  gas  flow  and  primarily 
broadened  Doppler  lines,  the  specific  expression 
of  nonsaturation  spectral  line-shapes  are  obtained. 
The  frequency-shift  relationship  and  the  qualitative 
and  quantitative  relationships  among  shifts  and 
deformations  of  nonsaturation  spectral  line- 
shape,  and  gas  macroscopic  flow-velocity  and 
macroscopic  flow-pattern  are  also  given. 


I.  Preface 


The  interaction  phenomenon  of  gas  flow  and  laser  is  somewhat 
different  from  the  interaction  phenomenon  of  gas  and  laser.  The 
macroscopic  flow  of  gas  will  cause  the  absorption  (or  emission) 
spectral  lines  to  have  frequency-shift  due  to  the  Doppler  effect, 
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and  deformation  of  nonsaturation  spectral  line-shapes  will  occur 

in  addition  to  frequency-shift  (saturation  spectral  line-shapes 

f  o  1 

will  certainly  also  have  a.  similar  situation  occurring  ) .  Frequenc 

shift,  spectral  line-shape  shift  and  change  in  line  width  after 

deformation  are  all  proportional  to  the  flow-velocity  component 

in  the  light-distance  direction,  thus  flow-velocity  can  be  measured 

according  to  laser  spectral  lines  measurements^^.  Parameters 

such  as  gas  flow  density,  temperature  and  gas  composition  naturally 

can  also  be  obtained  from  laser  spectral  lines  measurements. 

Using  spectral  lines  frequency-shift  to  measure  flow-velocity 

has  the  advantage  of  non-cc atact  measurement;  comparing  with  the 

*  [31 

laser  Doppler  velocity  measurement  method  (LDV)  ,  it  has  the 
advantage  of  not  requiring  the  addition  of  tracing  particles  in 
gas  flow;  therefore,  there  are  no  such  problems  as  disruption 
of  gas  flow  and  whether  particles  can  move  with  gas  flow,  nor 
are  there  limits  on  gas  flow-velocity  to  be  measured  (or  gas  flow 
Mach  number  Ma )  and  gas  flow  acceleration.  Spectral  lines  frequency- 
shift  is  proportional  to  the  flow-velocity  component  in  a  light- 
distance  direction.  The  larger  the  flow-velocity  is,  the  more 
effective  the  said  method  is.  Therefore,  it  has  the  unique  advantage 
to  transonic,  supersonic  flows  and  the  measurements  of  their  boundary 
layer  flows.  The  shift  and  deformation  situations  of  nonsaturation 
spectral  line-shape  are  related  to  the  gas  macroscopic  flow-pattern; 
therefore,  said  method  is  also  possible  to  be  used  i"n  diagnosing 
gas  macroscopic  flow-pattern. 

i  s  t  i  n  ;  u  i  ••  :i  ability  of  spectra  of  different  gas 
aoes  and  p  articles ,  as  well  as  the  development 


in  high  resolution  laser  spectral  technology,  it  is  possible  to 
more  effectively  utilize  the  laser  spectral  diagnosis  method  for 
parameter  measurements  of  mul t i -compos i t i on  gas  mixture,  isotope 
gas  mixture  and  multi-phase  flow,  and  meanwhile  obtain  parameter 
information  such  as  flow-velocity,  composition,  density  and  temperatur 
etc.;  .said  method  can  also  be  used  in  parameter  measurements  of 
various  inhomogeneous  gas  flow  media;  for  turbulent  flow  diagnosis, 
said  method  possesses  a  very  large  potential.  Hence,  it  appears 
quite  necessary  to  further  study  the  principle  of  spectral  velocity 
measurements,  to  conduct  analyses  and  discussions  on  the  situation 
of  spectral  lines  frequency-shift  as  well  as  shifts  and  deformations 

i 

in  spectral  line-shape. 

II.  Xonsaturat ion  Spectral  Line-Shape  of  Gas  Flow  Medium 


The  spectral  line-shape  definition  of  spectrum  is 

,  >>  ,  KO 

g(v,  v,)  —  —y-  (  2  .  •  ) 

/  —  j  l(v)dv,  j  ~g(v,  v,)dv  —  1 

where  v  is  frequency  is  spectral  line  center  freauencv;  7(v) 

J 

represents  the  distribution  function  of  light  intensity  based 
on  frequency,  I  is  total  light  intensity.  For  the  common  gas 
situation  in  fluid  dynamics,  the  contributions  of  Doppler  broadenin'; 
and  pressure  (homogeneous)  broadening  to  spectral  li'nes  broadening 
are  equally  important.  The  nonsaturation  spectral  1  me-  :h  e:e  - 
f  -j  r  urn  t-'t  o  tonic  flu  i  J  d  vnuni  c 


■;  \  a 


where 


,‘>(v>') "  sb  (“T  “»{-( "•'(gcr*)’!} 

Avh,I  2x 


Av0 

*»(*,  »')  “ 


f  2  .  5  > 


gD(v’,  vq)  is  the  Gauss  line-shapes  at  center  frequency  of  v  , 
g„(v,  v’)  is  the  Lorentz  line-shape  at  center  frequency  of  v1; 

A  v  and  A  v  are  the  f ul :idths  of  Gauss  and  Lorentz  spectral 
line-shapes  at  half-peal;  value,  respectively,  and  are  called  lino- 
shape  width.  At  low  gas  pressure  and  prcdomi na ntlv  Doppler  broadening 
(i.e.,  A  vD  » A  Vjj),  take  the  approximation  g^(v',  vQ)=gp(v,  vq  )  ,  and 
(2.3)  can  be  simplified  as  « 

?(k.  V,)  —  g[>Cv,  Vj)  (2.6) 

At  very  high  gas  pressure  and  primarily  homogeneous  pressure  broadeni 
(  i.e.  ,  A  vH  v  )  ,  take  the  approximation  g  (v,  v')  eg^fv,  v  ) 

and  (2.3)  can  be  simplified  as: 

g(v,  V,)  —  Zh(v,  V,) 


On  the  other  hand,  the  nonsaturation  and  absorption  (or  emission 
spectral  line-shape  of  gas  flow  can  be  theoretically  determined 
through  a  certain  velocity  distribution  function  of  particles. 

This  kind  of  particle  velocity  distribution  function  F(v’)  gives 
the  probability  for  the  particle  to  possess  a  displayed  velocity 
between  the  range  of  v'  and  v’-dV  along  the  light-distance  direction 
for  various  homogeneous  gas  flow  media,  the  displayed  velocity 
v ’  can  be  expressed  is  the  linear  summation  of  the  projection 
v  1 1  • ;  s  • : .  '  i  - 1  ,  2,  Vi  on  tn.e  1  i  air.  -  a  i  stance  direction  of  the  th  r  to 


coordinate  components  of  particle  thermal  flow-velocity  v„,  and 


local  average  gas  flow-velocity,  i.e., 

3 

‘  V  —  vT  4-  u*  (2.3  > 

•  »  I 

The  particle  velocity  distribution  function  F(v')  can  be  expressed 
as 

FW)  —  /,(«i)/r(vr)  (2.9) 

*-i 

where  f_,(v-p)  is  the  particle  velocity  distribution  function  correspon 
ing  to  particle  thermal  flow;  f  ^ ) ( i  =  1 ,  2,  3)  is  the  particle 
velocity  distribution  function  corresponding  to  gas  macroscopic 
flow; 


represents  the  combined  probability  for  particles  to  be  simultaneous! 
in  the  thermal  flow-velocity  range  between  v  and  v^-dv^,  and 
*  he  projection  values  range  between  u.  and  u^-du^(i  =  l,  2,  3)  on 
the  L i ght - a i stance  direction  of  the  coordinate  components  of  gas 
macroscopic  flow-velocity .  The  combined  probability  F  '  v '  )  is 
< "  ;  r  1 1 1  y  the  probability  for  the  particle  to  possess  a  displayed 
' 1  -  " : a  y  between  the  range  of  v ’  and  v'-dv’.  These  velocity  distribu 
■  .  functions  are  defined  as  the  unity  integral,  i.e., 

[  1r(vT)dvT  —  l ,  j  “  |  t,(u,)du,  —  1  '  _  .  .  O  ) 

n  -jjj  n  - 1  •  ' 2  •  - 1 } 

jjjj  1t(vt)  JJ  [1,(u,)du,]dv  T  —  1  (  2  .  '  2  '■ 

"  rJ” 


F(v")dv  —  1 


Here  u.  and  u.  are  the  two  boundaries  of  u.(i=l,  2,  3).  In 

10  1C  1 

order  to  solve  the  probability  function  F(v'),  the  following  mathe 

matical  transformation  is  performed  now 


Ut-VXu.,  vt)  (/  -  1,2,  3) 
s 

»  —  ~yi  «,  + 


(2.14) 


where  the  specific  function  form  of  LL(u^,  vT)  is  related  to  the 
specific  form  of  gas  macroscopic  flow,  but  the  following  equation 
is  also  true 


1  “  \\\\  /r(pr)  n  i /.(«.) du,] dvT 

-  jj]{  It «/..«0  n  UM, 


(2.15) 


where 
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9(1//,  O 
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dvT 
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ai/2 
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6Uj 

0H, 

d«; 

du, 

dvT 

at', 

au, 

dU, 

dU , 
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a«] 

Bu, 

dvr 

dr' 

at-' 

ar' 

Qv' 

(2.16) 


For  the  case  of  a  given  specific  flow  field,  the  integral  of  (2.14 
in  the  phase  space  ( Lh  ,  v')  can  be  converted  to  the  following 


form 


Therefore, 


«</,.,■)  n  >wi 

—  j  «>(Y,  «,»,  u.^dv'  —  * 


F(f')  —  <D(v’,  «,,) 


(2.17) 


(2.13) 


The  displayed  velocity  v'  of  particle  and  its  displayed  frequency 
v1  satisfy  the  Doppler  re  la i tonshi p ,  and  under  first-order  approxima¬ 


tion  of  — —  (<<  1  )  , 
c 


the  Doppler  relationship  is 

(»>'  —  Vo) 

V, 


(2.19) 


Thus  , 


(2.20) 


where  v  is  the  resonant  raidation  frequency  of  particle  at  rest, 
c  is  the  speed  of  light.  It  can  be  observed  from  ( 2 . 1 8 ) - ( 2 . 2 0 ) 


that  the  probability  for  the  particle  to  be  in  its  displayed  frequency 


range  from  v1  to  v'  *  dv '  in  the  frequency  space  is  also  the  frequency 

I 

distribution  function  of  the  particle,  which  is  — (v'-v  , 

o 

uio'  Uic^  and  satisfies  the  unity  integral  condition 


—  <P(u'  —  I',,  u,o,  <*,,)dv'  —  1 

v. 


For  the  general  situation  of  simultaneous  pressure  broadening 
and  Doppler  broadening,  the  computation  of  nonsaturation  absorption 
(or  emission)  spectral  line-shapes  must  not  only  compute  the  contribu¬ 
tion  of  the  particle  with  radiation  phase  resonance  at  frequency 
v  to  the  radiation  at  frequency  v,  but  must  also  compute  the  contribu¬ 
tion  of  the  particle  with  displayed  frequency  v'(v'^v)  to  radiation 
at  frequency  v  due  to  induced  homogeneous  broadening  effect. 

Therefore,  the  general  expression  of  nonsaturation  s’pectrul  line- 
shape  is 

gl(v,v,)—\  —  <P(v'  —  Vi,  ul!)g„{v' ,  v)dv' 

.  *3  U, 

For  the  soectral  line-shape  :i  ( v ,  v  '  of  snoot  r  a  1 

-  l  -‘  n 


corresponding  to  gas  unmacroscopic  flow  (refer  to  (2.3)),  there  are 
generally  frequency-shift  and  deformation,  broadening  or  narrowing 


of  line  width,  decrease  or  increase  in  peak  value.  Since 


A»d  -  ^  a,(  Ini)''1, 

c 


(2.23) 


where  a^  is  the  most  probable  thermal  f low-veloci ty ,  k  is  Boltzmann 

constant,  m  is  mass  of  particle.  When  gas  flow-velocity  is  of 

the  magnitude  of  the  speed  of  sound,  the  full  width  at  the 

half-peak  value  of  the  frequency  distribution  function  <^>(v'-vo, 

u.  ,  Uic^  '^S  also  t^ie  magnitude  of  doppler  line  width  ^VD- 

Although  the  displayed  velocity  v'  of  the  particle  might  be  different 

from  the  thermal  flow-velocity* vT  of  the  particle  (refer  to  (2.14)), 

3 

the  Z  u.  contained  in  v'  for  gas  microgroup  is  the  same  for 
i  =  1  1 

all  particles.  Therefore,  the  following  transformation  can  be 
performed 


v '  —  i*a  “  »'* 


—  v0, 


(2.24) 


where  7\  is  laser  wave  length,  vQf.  is  the  line-shape  center  frequency 
of  nonsaturation  spectral  line-shapes  of  gas  flow.  (2.24)  gives 
the  shifts  in  location  of  spectral  line-shapes  peak  value  caused 
by  gas  macroscopic  flow.  The  frequency-shift  (v  -v  )  is  positive 

J 

or  negative  based  on  whether  the  direction  of  u-  is  concurrent 

•  •l 

or  countercurrent  to  the  laser  propagation  direction.  The  nonsatura¬ 
tion  spectral  line-shapes  of  gas  flow  is  symmetrical  or  near- 
symmetrical  to  v=v  .  Considering  (2.24)  and  also  referring  to 
the  following  text,  (2.22)  can  be  transformed  to 

£/(«*>*'«)“*(  —  <P(v",  v)dv"  (2.25) 


WO  • 


low  or' 


;  u  r  c 


a  n  c! 


;  a  s 


low  spectral 


l  ne 


shapes  width  A  >>AvH,  the  approximation  of  ^  (v'-vo,  u^q, 

u.  )=  3)  (v-v  ,  u.  ,  u.  )  and  u>(v",  v  J.)=${v,  v  ,)  can  be  adopted 
ic  o  10  ic  2  of  X  of 

in  performing  the  integral  of  (2.22)  and  (2.25).  Therefore,  from 
(2.22  )  and  (2.25) 


g/(v,  W0)  “ 1  <£(w  Wo,  M,o»Mi») 

Wo 

£,(w,  Wo/)”  —  /) 

Wo/ 


(2.26) 


(2.27) 


Generally  speaking,  —  ^(w,  r0/)  ^c(w,  k,/).  For  the  case  of  very 

w„ 

high  gas  flow  pressure  and  AvH»A^  ,  the  approximation  of 
g„(v',  v)=g„(v  ,  v)  and  g„(v",  v)=gtI(v  , ,  v)  can  be  adODted  in 

n  HO  H  H  O  r 

performing  integral  of  (2.22)  and  (2.25).  So,  from  (2.22)  and 
(2.25) 


f /(*»  «•)  j'») 


(2.28) 


i/(w»  »«)  ~  *«(»»♦  *ty) 


(2.29) 


(2.28)  and  (2.29)  are  consistent  with  (2.27),  but  they  also  contain 
a  layer  of  new  significance,  i.e.,  under  high  gas  flow  pressure 
and  Av^  »  A  ,  the  nonsaturation  spectral  line-shapes  of  gas 
flow  is  independent  of  gas  flow  conditions,  and  they  are  all 
homogeneous  Lorentz  line-shapes. 


In  addition,  under  the  condition  of  u^(i=l,  2,  3)  approaching 
0,  i.e.,  gas  macroscopic  flow,  it  is  known  from  (2.14)  and  (2.15) 
that  the  particle  frequency  distribution  function  approaches  Gauss 

line-shaoes,  i.e.,  —  9{v'  —  w„  0,0)  -**»(*’,  w,),  ancj  therefore,  the  nonsatura 

‘  »■«  ' 

tion  spectral  1  ine-sahpes  equation  (2.22)  is  simplified  to  spectral 


III.  Derivation  and  Example 

We  will  conduct  specific  discurssions  below  on  the  three 
typical  gas  flows  and  the  situation  with  predominently  nonhomogeneous 
Doppler  broadening 

1.  Shear  Flow  and  One-Dimensional  Homogeneous  Flow 

Assume  the  whole  body  (i.e.,  macroscopic)  of  a  gas  microgroup 
which  is  small  macroscop i cal ly  and  large  microscopically  near 
a  certain  point  in  space  moves  along  the  x  direction,  and  the 
magnitude  of  velocity  is  propottional  to  y  (refer  to  Fig.  1),  i.e., 

v,  -  y  +  u„v,  —  I/,  -  o  (3.1) 

where  1  is  the  characteristic  line  degree  of  the  gas  microgroup. 

(3.1)  also  represents  the  flow  along  the  x  direction  of  the  entire 

gas  flow.  The  velocity  gradient  along  the  y  direction  is  a  constant 

r  5  1 

shear  flow,  e.g.,  Couette  flow  ,  etc.  The  special  case  of  U  =U  , 

G  O 

U  =U  =constant  in  the  shear  flow  equation  (3.1)  is  a  one-dimensional 
homogeneous  flow.  Assume  the  laser  measurement  is  conducted  along 
A  and  B  directions  (refer  to  Fig.  1);  A  direction  is  perpendicular 
to  gas  flow  direction  (i.e.,  x  direction),  and  B  direction  is 
at  an  angle  9  with  A  direction.  For  B  direction,  the  projection 
values  of  the  three  coordinate  components  of  gas  flow-velocity 
are,  respectively 


Fig .  1 .  Shear  flow 

Key:  (1)  A  direction;  (2)  B  direction- 


considering  the  fact  that  gas  flow  density  is  constant  and 

velocity  distribution  is  linear  distribution;  therefore,  the  velocity 

distribution  function  f.(u.)  is 

l  l 


/,(«.) 


1(m,  —  U,  Sin'S),  V,  -  Ut  -  U. 
‘  1 

(l/, -U,)s in  (3  U‘>U> 

/,(«,)  -  B( 0),  Uus)  -  B( 0) 


(3.3) 


^  is  a  function  of  £  .  The  following  mathematical  transformation 
is  performed 

L\  “  U-t  Ut  —  *1,  U ,  —  v'  —  tlj  •+•  U,  /  -5  -  \ 


7hen 


whe  re 


£(«■,  gr)  _  _J_ 

0(i/„  *0  2 

jj(j  /Krr)  n 

f.  - 1 - 

J--  2  (l/,  —  l/,)  sin  9  l  V  a,  J 

2  r  * 

:  —=  |  :  ,  the  probability  intoqn: 


(  3  .  o  ) 


l  .5 


assumed  to  be  the  Maxwell  balanced  distribution  of  a  5  mule  voioci 
component  when  (3.5)  is  derived.  Using  (5.19)  and  (5.50',  the 
nonsaturation  spectral  line-shapes  of  shear  flow  are  obtained 
from  (3.5) 


„)  -  -  [erf  (lV^~2 

*'  ’  Av0  2 ([/,-£/,)  sin  d  t  V  Avo  <,  / 

-  -rf  (ly/^i  iUZLfi  +  ‘  ' 

\  Art,  f  a,  /l 

When  9  =  0,  i  .  e  .  ,  investigating  along  A  direction,  the  nons.it  ur  at  io 


spectral  line-shape  is  regressed  from  (3.6)  to 


(  3.7) 


This  is  the  Gauss  line-shape  that  coincides  with  the  nonsaturation 

spectral  line-sahpes  of  gas  macroscopic  flow.  For  one-dimensional 

homogeneous  flow,  i.e.,  when  U  =const®ht  and  U  -’J  ,  (3.6)  is  re:; re 

x  e  o 


*/("»  »it) 


2  la2 

A  v0  V  * 


—  v,  —  Ai Avt  —  —  -  Av°  sin.0«>O.6M  j,Avd  >in9 

•>  2  Via  2 


(3.3) 


Aw.  —  M-fs  jin  9 
1 


(3.9) 


•.’here  Ma  is  the  gas  flow  Mach  number.  (3.3)  and  (3.°)  show  that: 
for  one-dimensional  homogeneous  :  s  flow  ( U  __  =  U  =U_. )  ,  the  nonsatura 
spectral  line-shapes  measured  along  3  direction  are  still  Gauss 
line-sahpes,  and  said  line-shapes  have  f  requency  -  sh  i.f  f  with  re  spec 
to  the  Gauss  line-shapes  measured  along  A  direction.  The  amount 
of  frequency-shift  A  v  of  the  location  of  1 ino-<nhpos  coax  v  1T1  1 
is  given  by  (3.9).  These  conclusions  are  cons :  mvi'  ‘  h 

analytical  conclusions  of  7e Terence  5  '.  >  :t  1  >  f  • 


[2]  indicate:  if  the  time  required  for  the  gas  flow  to  pass  through 

a  laser  light-beam  cross  section  can  be  compared  with  molecular 
energy-level  collision  life  ( e .  :j .  ,  when  light-beam  width  is  0.:  cm 
and  flow-velocity  perpendicular  to  light-beam  is  ax'.  0~  m/sec, 
the  time  for  gas  flow  to  pass  through  the  light-beam  cross  section 
is  only  2x10  ^  sec),  then  the  gain  saturation  rule  of  gas  flow 
medium  and  saturation  spectral  line-shapes  will  all  be  different 
from  the  situation  of  gas  unmac roscop i c  flow.  However,  the  non¬ 
saturation  gain  of  gas  flow  medium  and  nonsaturation  spectral 
line-shapes  are  consistent  with  the  situation  of  gas  unmacroscopic 
flow;  the  aforementioned  analyses  reconfirm  this  conclusion  and 
further  show  that:  the  nonsaturation  spectral  line-shapes  measured 
leaning  toward  the  gas  flow  direction  can  be  obtained  through 
frequency-shift  changes  of  nonsaturation  spectral  line-shapes 
measured  perpendicular  to  the  gas  flow  direction.  The  amount 
of  f requencv-shi ■ t  A  v^  of  the  location  of  spectral  line-shapes 
peak  value  is  given  by  (3.9).  A  vf  is  proportional  to  the  product 
of  f low-ve loc i t y ,  A v^  and  the  sine  of  the  angle  between  the  two 
directions,  and  Fig.  2  gives  the  relationship  of  Av  /A  v„  with 
the  change  of  gas  flow  Mach  number  Ma  ;  A  vf  is  inversely  proportion) 
to  wave  length  ^  ;  refer  to  equation  (3.9).  This  is  consistent 
with  the  relationship  of  single  spectral  line  Doppler  frequency- 
shift.  Therefore,  in  the  domain  from  the  infrared  to  ultraviolet 
laser  spectrum,  a  gas  flow-velocity  of  1  m/sec  will  cause  a  frequency 
shift  between  3 .  1  Ml!  and  '.0  MH  approximately. 
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Fig.  3.  Two- dime  ns 
expansion  flow 
Key  :  (  1  )  A  direct: 

(2)  3  direction. 


For  shear  flow  condition  of  U  equation  (3.6)  can  be 

t 

,  r ins  formed  to 


L'f  ~  U . 

2a, 


2s/ In  2  fl> 

Av0  2(f/, —  U,)an8 


erf  I  1  \/  lc  2 


;  »  —  Oo  —  f/) 


sindj  —  erf  ^2\/ In  2 


AvD 

v  ~  ( t.0  —  )  i/,  —  L  „  . 


Al-J 


-a. 


w  here 

_  U,  4-  V,  _Ar£_  (  3  .  ;  i  > 

l;  :-v/ln2 

If  can  be  observed  from  (3.10)  and  the  following  text  that  she 
flow  nonsaturation  spectral  line-shapes  are  equaivalent  to  non 
f ion  spectral  line-shapes  of  two-dimensional  symmetrical  exp  an 
f  1  o  w  . 


2.  Two-Dimensional  F.xpansion  Flow 


In  the  domain  of  detection  points,  as  sumo  t  f ;w-vol 
component  along  axi  a  1  direction  '  d  i  re  •-  i  on  •  -  •  x  • 

/elocify  component  •  i  lor:  g  1  i  *■  e  r  a  1  i  i  re-**-  •;  > :  •. 


linearly  from  -w  to  wg ;  refer  to  Fig.  3 ;  for  the  two  detection 
directions,  A  direction  is  perpendicular  to  x-axis,  3  direction 
is  at  an  angle  0  with  A  direction.  The  projection  values  of  the 
three  coordinate  components  of  flow-velocity  are 


*i  “  V ,  sin  6  *«  V c  sin  8,  — »  9,  us  “  u  cosd 

w,  —  w, 

w  “  — - — 2  z  -4-  uj. 


(3.12) 


For  homogeneous  expansion  flow,  it  can  be  assumed  that  within 
the  velocity  range  of  -v  $  W  $  Wg  the  distribution  probability 
of  gas  particles  in  the  macroscopic  flow-velocity  space  is  constant 
so  the  velocity  distribution  functions  of  particles  are 

t 

—  5(«,  —  t/,sin 6),  /;(«,)  —  5(0) 

5(°)  — o 

~  +  ^  g  l«,|  <  max(  |  If'd  cos6,  |W’,]cos0) 

0  |«il  >  max(  |  H't  j  cosd,  |H',Jcos3) 


Perform  transformation 


L‘  -  U,  —  u„  U , 


f'r  .v.  i  «  t  »,  t 


llj/rC^n  U,(u,)au,  ]ct - 


(3.14) 
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!  —  2(  tv3  4  W,)cosd  l"  y  ap  ) 
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where  v0f=v0~Avf'  see  equation  (3.9).  In  equation  (3.16),  let 

A  vf=0  and  obtain  the  nonsaturation  spectral  line-shapes  measured 

along  A  direction.  So,  the  constant  velocity  macroscopic  flow 

of  gas  along  x  direction  will  cause  the-  nonsaturation  spectral 

line-shapes  measured  in  B  direction  to  have  a  specific  frequency- 

shift  with  respect  to  the  spectral  line-shapes  in  A  direction. 

The  amount  of  f requ f ency-shi f t  is  consistent  with  the  amount  of 

line-sahpes  center  frequency-shift  A  v^  of  gas  unmacroscopic  expansion 

flow  (see  (3.9)),  and  this  indicates  that  gas  with  nonsaturation 

spectral  line-shapes  of  complex  macroscopic  flow  is  the  superposition 

of  gas  with  spectral  line-shapes  of  simple  macroscopic  flow. 

If  gas  flow  expansion  is  weak,  W  /a  «1,  W  /a  <\  1  ,  then  equation 

e  p  op 

(3.16)  is  regressed  to  equation  (3.8);  if  the  effects  of  gas  flow 
macroscopic  expansion  flow  cannot  be  neglected,  the  expansion 
flow  will  make  nonsaturation  spectral  line-shapes  become  non-Gauss 
line -shapes,  and  the  spectral  line -shape  peak  value  decreases 
(expansion  condition)  or  increases  (convergent  gas  flow  condition), 
line  width  broadens  or  narrows.  The  location  of  the  spectral 
line-shape  peak  value  obtained  from  (3.16)  is 

y.  -  y»  -  Ay,  +  -  W ,)  (3.17) 

w  4Vln2  "»■ 

Therefore,  the  amount  of  frequency-shift  Avg  of  asymmetrical 
expansion  spectral  line-shape  peak  value  location  with  respect 
to  symmetrical  expansion  <  K  -W  >  spectral  line-shape  peak  value 


4  V  In  2  a  t 


Substitute  (3.18)  into  (3.16)  to  obtain  the  magnitude  of  spectra: 
line-shape  peak  value 


OteXgl(v,  V')  —  g,(  vm  9  *j) 

2-vAn? 


r  cos  9  i 


The  half  width  A  v*  =  ( vD-vof. )  at  the  half  peak  value  of  spectral 

line-shape  is  obtained  from  the  relationship  below;  is  the 

frequency  at  half  oeak  value  location,  v=v  -  v*, 

^  *■  D  of  D 


erf ‘f 2 V la  2  ^S.  +  *!«**)  '-1  er4,  ( 2sJ\a2  -  *£S5«l) 
\  Avp  /  *  \  Ai>d  / 

j  i  7  • 

-  erf  [■£!£(  +  ^)] 


(3.20) 


T,\'e  better  make  some  generalized  mathematical  discussions  on  expansion 


flow  nonsaturation  spectral  line-shapes,  let 

Wt  -  aW c,  X  -  2  >/to2  v-^t,  Y  -  i^25» 

Ai;0  a. 


,  x  Avd  v  * 

h.  ■ 


_-erfU 


(3.21) 


(3.22) 


Perform  mathematical  transformation 


*'  +  Y’  -  x  +  y,  x'  -  r  -  x  -dv 


3.23) 


e(X,  Y)  -£(X\  ?’)  -  ^-Ierf(X'  +  V')  -  erf(X'  - 
■  -  -  ■ '  $  j?  —  X’  —  y  _  — h —  Y' 

'  .  1  +  •  1  •+•  0  . 


(3.24) 


m*xKX,  y)  -  nu xg(X',  r)  -  -^TTT  erf(V') 


■■  3  .  3  “  i 


flow  can  bo  obtained  from  the  spectral  line-shapes  of  expansion 

flow  through  the  near-jet  transformation  equation  (3.23).  From 

(3.10)  it  is  also  known  that  shear  flow  spectral  line-shapes  are 

also  equally  equivalent  to  the  spectral  line-shapes  of  a  symmetrical 

expansion  flow.  Figure  A  gives  several  spectral  line-shapes, 

where  curve  1  is  a  Gauss  line-shape;  curve  2  is  a  symmetrical 

expansion  flow  (K  =  W  )  spectral  line-shaoe  and  curves  3  and  -4 

o  e 

are  asymmetrical  expansion  flow  spectral  line-shapes.  Figure  3 

gives  the  relationship  curves  of  nonsaturation  spectral  line-shapes 

oeak  values  with  the  changes  in  W  cos9/a  .  Peak  values  droo  monot- 

J  e  p 

onouslv  as  Wecos/a  increases. (  See  Fig.  6  for  the  relationship 
of  line  width  nonsaturation  spectral  line-shapes  with  the  changes 
in  V.'ecos9/ao .  The  increase  in  line  width  is  comparable  to  Av^; 
therefore,  the  expansion  flow-pattern  and  magnitude  of  expansion 
velocity  can  be  measured  according  to  the  changes  in  line-shapes 
a  n  d  line  width. 

Tn  the  relationship  equation  (3.15)  of  expansion  flow  spectral 

1 i.oe-shaoes ,  if  we  let  W  =W  and  use  (U  -U  )  s  inQ/2  to  substitute 

e  o  e  o 

V.'  cos9,  then  equation  (3.16)  is  transformed  into  shear  flow  spectral 
line-shapes  equation  (3.10).  Hence,  the  analytical  conclusions 
md  those  relational  curves  of  symmetrical  expansion  flow  in  Fig.  4 
through  Fig.  6  are  equally  suitable  for  shear  flow.  Here  we  omphasi 
again  that:  for  emission  spectral  1  ine-sahpes ,  the  flow  indicated 
in  the  above  analyses  is  the  flow  in  the  neighboring  r on  ion  of 
the  riaiation  emission  point  (i.e.,  measurement  point;,-  therefore , 


Closing  Remarks 


Spectral  lines  and  nonsaturation  spectral  line-shpaes  will 
have  frequency-shift,  deformation  and  change  in  line  width  due 
to  macroscopic  flow  of  gas;  within  the  domain  from  the  infrared 
to  ultraviolet  laser  spectrum,  a  gas  macroscopic  flow-velocity 
of  1  m/sec  will  cause  f requency-shi f t  and  line  width  change  between 
0.1  MH  and  10  MH  ,  and  these  frequence-shifts  and  changes,  in 
turn,  can  be  used  to  measure  macroscopic  flow-velocity  and  flow- 
pattern  of  gas.  This  measurement  method  has  advantages  such  as 
noncontact  measurement:  no  need  to  add  tracing  particles  which 
are  far  larger  than  molecular  volume  in  the  gas  flow;  no  interference 
with  flow  field;  no  limitations  to  magnitude  of  flow-velocity 
and  flow-acceleration;  being  able  to  conduct  simultaneous  measurement 
of  multiple  points;  and  being  able  to  obtain  parameters,  such 
as  gas  density,  temperature,  compositions,  etc.  simultaneously. 

It  also  possesses  great  potential  for  measuring  various  nonhomogeneous 
gas  flows  and  turbulent  flows. 
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